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Recently, the quantum spin Hall effect has been predicted in (11 l)-oriented thin films of SnSe and SnTe topo¬ 
logical crystalline insulators [1,2]. It was shown that in these films the energy gaps in the two-dimensional band 
spectrum depend in an oscillatory fashion on the layer thickness - the calculated topological invariant indexes 
and edge state spin polarizations show that for films 20-40 monolayers thick a two-dimensional topological 
insulator phase appears. Edge states with Dirac cones with opposite spin polarization in their two branches 
are obtained for both materials. However, for all but the (11 l)-oriented SnTe films with an even number of 
monolayers an overlapping of bands in T and M diminishes the final band gap and the edge states appear either 
against the background of the bands or within a very small energy gap. Here we show that this problem in SnSe 
films can be removed by applying an appropriate strain. This should enable observation of the Quantum Spin 
Hall effect also in SnSe layers. 


I. Introduction 

Quantum spin Hall (QSH) effect is a net result of oppositely 
polarized spin currents traveling in opposite directions along 
the edges of a two-dimensional topological insulator (2D TI). 
In 2D TI the metallic spin-polarized edge states are protected 
by time-reversal symmetry and robust against back scattering. 
Therefore the spin-polarized currents can flow without dissi¬ 
pation. This is an extremely appealing feature for low-power- 
consumption usage, which triggered an ever growing effort to 
establish system that harbors QSH phase. 

Lately, it has been theoretically predicted that for (111)- 
oriented thin films of SnTe and SnSe of particular thicknesses 
the QSH state can be obtained [1,2]. Both these materials are 
topological crystalline insulators (TCIs) whose surface states 
are protected not by time-reversal but by mirror plane symme¬ 
try [3, 4]. When the surfaces of the layer are well-separated, 
they exhibit four Dirac cones at the four projections of L 
points onto the 2D Brillouin zone of the (111) surface (one 
in T and three in M points) [5, 6]. Otherwise, in ultrathin 
films, the surface states of top and bottom surfaces start to in¬ 
teract and energy gaps open where the Dirac nodes used to 
be. These energy gaps can be modulated by changing the film 
thickness - it was shown that for some thicknesses the 2D va¬ 
lence and conduction bands intersect and the 2D states energy 
gap reaches negative values. In SnTe films only the one energy 
gap at T changes sign for some thicknesses, while in SnSe the 
three gaps at M points oscillate with thickness. Thus, in both 
SnTe and SnSe layers the band inversions take place at odd 
number of points. In Ref. [1] it was shown that indeed a topo¬ 
logical phase transition from trivial insulators to 2D TIs can 
be achieved by tuning the layer thickness. It is worth mention¬ 
ing that a similar, thickness-induced phase transition was ob¬ 
tained also for (OOl)-oriented ultrathin film, however, in this 
case the conversion is between trivial insulator and 2D TCI 
phase [7, 8], 

As it was discussed in Ref. [ 1 ], for all but the (111 )-oriented 
SnTe films with an even (ca 20) number of monolayers, the 


QSH effect would be difficult to observe because the T point 
and one of the M points are projected to the same k x = 0 point 
in the ID BZ of the [110] edge. Due to the overlapping of 
these two projections the final band gap can be considerably 
diminished by mutual shadowing. As a result the edge states 
appear either against the background of the bands or within a 
very small energy gap. In this article, we scrutinize the possi¬ 
bility of reducing such band gap shadowing in (111 [-oriented 
SnSe films by biaxial strain. We predict that applying appro¬ 
priate strain, e.g., by well chosen substrate, should enable ob¬ 
servation of the QSH also in SnSe thin layers. 

II. Band structure in strained layers 

We investigate the effect of biaxial tensile strain on the en¬ 
ergy spectrum, in particular on energy gaps, of (11 l)-oriented 
SnSe film. The film lays in the x-y plane and the z direction is 
chosen along [111] crystallographic axis. It is well known that 
strain changes the electronic properties of the material. It can 
shift and split the degeneracy of energy levels. In heterostruc¬ 
tures strain affects also the staggered band offsets and is con¬ 
sidered as a powerful tool to engineer alternative routes for 
enhancing semiconductor structures [9] . For a homogeneous 
semiconductor, the variation of energy gap as function of the 
imposed strain is qualified in terms of deformation potentials 
tensor (D) and strain component (e). When the deformation 
is not very strong and the distance between atoms is still close 
to the equilibrium lattice constant a 0 , D can be regarded as a 
constant coefficient and hence, the change of the energy with 
strain is approximately linear. Under isotropic biaxial strain in 
x-y plane, longitudinal s\\ and transverse e ± strain components 
are expressed as following: 

_ _ clx ’ y _ _ _ Q z ,,, 

ao ao cio 

where «d = «n( 1 + deformation^), a- = ao( 1 + deformation,.) 
and «o is the relaxed lattice constant. The longitudinal and 
transverse strain components are related through the Poisson’s 
ratio v that characterizes the elasticity of a crystal: 
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The Poisson’s ratio of rock-salt SnSe is assumed to be similar 
to the one for PbSe and taken from Ref. [10] to be equal to 
1.16 at low, T a 4 K, temperature. The lattice constant ciq — 
5.97A of the relaxed cubic SnSe was obtained in Ref.[ll] by 
first principle calculations. 
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FIG. 1. Calculated band structure for relaxed (red dotted lines) and 
1.14% deformed (solid blue lines) anion-terminated 17-monolayers 
of SnSe. 


In Fig. 1 we present the band structure of strained and fully 
relaxed 17-monolayer thick, anion-terminated (11 l)-oriented 
SnSe slab. The calculations were performed within a tight- 
binding approach. We have chosen this film to study the pos¬ 
sibility of aligning the energy gaps by strain, since here the 
overlapping practically closes the energy gap. Indeed, as one 
can see in the figure, in the relaxed layer the energy gap at 
M is completely covered by the valence band in the vicinity 
of T point. Our calculations show that while the deformation 
increases, the energy gap at F point is shifted down in energy, 
whereas the band gap at M point is moving upward, as shown 
in Fig. 2. This is due to the difference between the effective 
masses in the vicinity of these two points of the Brillouin zone 
and subsequently different deformation potentials. 

For deformation equal to ~ 1.14%, the tops of the valence 
band at I and M have the same energy (compare Fig. 2) and 
the two energy gaps are aligned. As a result the band gap in 
the 2D spectrum of such slab opens - the final total energy 
gap is equal to ~ 0.57 meV, as shown in Fig. 1 by the green 
strap. 

The lattice constant associated with the optimal defor¬ 
mation, which aligns the energy gaps in the anion-terminated 
17-monolayers of SnSe, is equal to 6.04 A. Thus, the epitaxial 
growth of SnSe layers on a substrate whose lattice constant is 
equal to 6.04 A would lead to the desired deformation. The 
appropriate substrates that can provide such lattice mismatch 
are for example PbS,Se|_ A with x-0.46 and PbS v Tei_ v with 
x^0.8. 



FIG. 2. The dependence of the energy gap of anion-terminated 17- 
monolayers of SnSe under biaxial strain on the percentage of defor¬ 
mation. The extrema of the valence band at T and M are at the same 
energy level when the structure is ~ 1.14% deformed. T v , r c , M v 
and Me stand for valence and conduction band energies at T and M 
points, respectively. 


III. Edge states and their spin polarization 

According to our previous results [1], a 17-monolayers 
anion-terminated SnSe film should be a non-trivial insulator 
whose band gap inversion occurs at three M points. There¬ 
fore, at the projections of the M points on the [110] edge, 
i.e., at k x = 0 and k x - ± VZr/a points of the ID BZ, three 
Dirac nodes are expected. However, the observation of these 
edge states would not be possible in the relaxed structures due 
to the mentioned above problem of strong shadowing of the 
band gap in such layer. Here, to endorse our strain-induced 
alignment solution, we have calculated the [110] edge spec¬ 
tral functions for both, relaxed and deformed, (11 l)-oriented 
anion-terminated SnSe films consisting of 17 monolayers. 
Fig. 3a, shows the edge states in the vicinity of k x = 0 for 
the relaxed structure. Notably the band gap is covered with 
the 2D valence band from F point. In turn. Fig. 3b shows the 
edge state of the same film, which is now 1.14% stretched. It 
evidently verifies that applying appropriate strain removes the 
shadowing of the gap. Now the edge state with a Dirac node 
can be clearly seen within the energy gap. 

Ultimately, we have calculated the spin-up (GT) and spin- 
down (GJ,) contributions in the edge spectral function of the 
17-monolayers anion-terminated SnSe film with 1.14% defor¬ 
mation. The spin polarization of the edge states is defined as 
G| - GJ,. The results are presented in Fig. 4. As expected 
all three Dirac crossings with opposite spin-polarization in the 
two branches are now located inside the band gap. This should 
enable their observation. 

IV. Conclusion 

In this article we show that biaxial strain can change the rel¬ 
ative positions of various energy gaps in the 2D band spectrum 
of (11 l)-oriented SnSe ultrathin layers. For all such films, an¬ 
ion and cation terminated as well as films with even number 
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FIG. 3. The calculated spectral function of (a) relaxed and (b) 1.14% 
deformed, 17-monolayers anion-terminated SnSe film at k x = 0. 
(The black solid lines denote the projection of the 2D bands bound¬ 
aries.) 
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FIG. 4. Spin density of 1.14% deformed 17-monolayers anion- 
terminated SnSe film in the vicinity of (a) k x = VT/t/c/d and (b) 
k x = 0. (Red and green lines denote the spin-down and spin-up po¬ 
larization, respectively.) 


of SnSe, cation-terminated, requires 1.07 % deformation to 
achieve an optimal alignment of the gaps, while 0.8 % defor¬ 
mation is needed to cause the same in a 20 -monolayer thick 
SnSe film. This means that the in-plane lattice constant of the 
film should be expanded to the value of about 6.02-6.04 A de¬ 
pending on the number of monolayers. The films can be thus 
stretched by using a substrate with such lattice constants. The 
substrate should be also a trivial insulator with a wide band 
gap to ensure that the interface between the SnSe layer and 
the substrate simulates as much as possible a free surface. It 
seems that the best candidate for the substrate is PbS,'fe|„ f 
with x-0.8 or PbS. t Sei_ v with x-0.46 for films with an odd 
number of layers and with x-0.56 for 20 -monolayer thick 
SnSe film. The energy gaps in the mentioned above com¬ 
pounds are ca 2-2.5 times bigger than in the film with odd 
number of monolayers, i.e., the band gap in the 17-monolayer 
thick SnSe slab is equal to -0.1 eV while for PbSo. 46 Seo .54 
£g=0.21 eV and E g = 0.27 eV in PbSo.sTeo. 2 . In the even num¬ 
ber of SnSe monolayers case the difference between energy 
gaps of the film and substrate is larger - for 20 -monolayers 
the gap is equal ca -0.04 eV, what means that the energy gap 
0.27 eV of PbSo. 56 Teo .46 is nearly 7 times bigger. It should 
be mentioned, however, that sulfur is usually not allowed in 
MBE systems and to grow such substrates other methods have 
to be applied. Finally, we note that the appropriate strain can 
be further adjusted by using the difference in thermal expan¬ 
sion coefficients of the SnSe and the used substrate. Of course 
one can think also about the possibility to obtain the needed 
deformation, i.e., the biaxial in-plane extension, by applying 
mechanically uniaxial pressure perpendicular to the ( 111 ) sur¬ 
face of the layer. 
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